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Abstract 

The A-dependence of the saturation momentum and the scaling behav- 
ior of the scattering of a small dipole on a nuclear target are studied in 
the McLerran-Venugopalan model, in fixed coupling BFKL dynamics and 
in running coupling BFKL dynamics. In each case, we find scaling not too 
far from the saturation boundary, although for fixed coupling evolution the 
scaling function for large A is not the same as for an elementary dipole. We 
find that is proportional to A^^^ in the McLerran-Venugopalan model 
and in fixed coupling evolution, however, we find an almost total lack of 
A-dependence in in the case of running couphng evolution. 

1 Introduction 

The study of high density QCD, states or systems where gluon occupation 
numbers are large compared to one, has become one of the central topics 
in QCD. There is a considerable literature whose focus is explaining much 
of small X and moderate HERA data[l-4] as well as the general features 
of hadron production in ion-ion collisions at RHIC[5-11] using high density 
wavefunctions as the basic ingredient. In any such description the central 
parameter is the saturation momentum, Qs, the scale at and below which 
occupation numbers large as 1/a in the light-cone wavefunction. 

The energy dependence of the saturation momentum has been widely 
studied[12-17] and it now appears that one has pretty good control over that 
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dependence in the context of rcsummed next-to-leading order BFKL dynam- 
ics. Recently there has also appeared some discussion of the A-dependence 
of the saturation momentum for large nuclei[ll-18] beyond the McLerran- 
Venugopalan model[19] model. The conclusion has been pretty much the 
same as for the McLerran-Venugopalan model, that is Ql{A) ~ A^''^. 

Our purpose here is to revisit the question of the A-depcndcncc of the 
saturation momentum implementing the BFKL dynamics more fully than has 
been done so far. In case fixed coupling BFKL dynamics is used we confirm 
that the primary A-dependence of Q^(^) is proportional to A^^^. The result 
is given in (26). We also confirm that the scattering of a small dipole, of size 
1/Q, on the nucleus gives a scattering amplitude which is only a function of 
Q'^/Q^A) so long as one is not too far from the saturation boundary. We 
note, however, that the scaling function is not the same as for the scattering 
of a dipole on another elementary dipole. (In the fixed coupling problem 
it does not appear possible to sensibly talk of dipole-proton scattering in a 
pcrturbativc context.) The results for the scattering amplitude are given in 
(27) and (29). 

In case a running coupling BFKL dynamics is used the situation changes 
radically[16,17]. Now the saturation momentum becomes almost completely 
independent of A at very large rapidities. Eq.(44) gives our result here for 
the weak A-dependence which is present. 

In both the fixed coupling and running coupling situations we have found 
it useful to view the BFKL-QCD evolution in an unusual way. The traditional 
way to view the scattering of a small dipole on a larger dipole, or on a hadron 
or nucleus, is to view evolution as part of the wavefunction of the larger object 
which is then probed by the smaller dipole. In the present problem it seems 
more convenient to view the BFKL evolution as part of the wavefunction of 
the smaller dipole, which evolution produces gluons (dipoles) of larger scale 
which eventually interact with the hadron or larger scale. 

Finally, because is not very large we do not feel confident in deciding 
which picture, McLcrran-Vcnugopalan, fixed coupling BFKL dynamics or 
running couplng BFKL dynamics is more appropriate for RHIC energies 
and for the LHC heavy ion regime. It might well be that the McLerran- 
Venugopalan model along with a modest amount of fixed coupling evolution 
is dominant at RHIC energies so that ~ A^^^ is the appropriate behavior 
of the saturation momentum. We would guess that the running coupling 
regime is likely dominant for protons and small nuclei at HERA energies. 
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and is perhaps important for large nuclei at LHC energies, but this is by no 
means certain. It would be very useful to have some numerical calculations 
which are reliable near the saturation region to try and see what the A- 
dependence of Qg actually is in the energy regions covered by the different 
accelerators. 



2 The semiclassical region (McLerran- 
Venugopalan model) 

We begin our discussion of the A-dcpcndence near the saturation boundary 
with a brief review of the McLerran-Venugopalan model[19]. Let T^r be the 
amplitude for the scattering of a quark-antiquark dipole of size Axj, = 1/Q 
on a nucleon. Using a normalization where the cross-section for scattering is 

On = 2rjv (1) 

one has 

^iv = (2) 

where xG^ is the gluon distribution of the nucleon evaluated at scale Q^. 
Let T^ip) be the corresponding amplitude for the scattering of a dipole of 
size on a nucleus at impact parameter h. Then, the dipole-nucleus cross- 
section is 

<^A = 2j d%TA{b) (3) 

with[19,20] 

T^^l-expi-^] (4) 



where 



QUA) = 2aM ■ 2y/W^ pQ' = jfQ'{A) (5) 



with, as usual[21] 
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Ql (A) = 2VW^^pxG. (6) 



c 



p is the nuclear density, Qg is the quark saturation momentum and Qs the 
gluon saturation momentum. It should be emphasized that Q = defines 
Q. Q is not a true momentum variable connected to x^— dependence by a 
Fourier transform. 

Eq.(4) has all the basic properties which we wish to note in this section. 
First of all (4) exhibits geometric scaling[3]. That is Ta is only a function 
of Q'^g{A)/Q'^.TA exhibits saturation so that ~ 1 when Qy{4Q'^) > 1. 
Finally, when 4Q^/Ql ^ 1,Ta is additive in the various nucleons. That is 

CTiV 



Ta ^ '^^2VW^p, (7) 
SO that there are no non-trivial nuclear effects present when Q"^ IQ\ ^ 1. 



3 Fixed coupling BFKL dynamics 

Now we shall consider the scattering of a dipole of size Ax^ = 1/Q on a dipole 
of size Aa;_L = (where Q'^ j p? ^ 1) and on a "nucleus" made up of A 
such dipoles distributed uniformly in a sphere of radius R and having density 
p — AI{^B?\ Of course this is not a realistic nucleus, since it is not made 
up of realistic nucleons, but it is not possible to deal with a realistic nucleon 
in the fixed coupling regime. Since our object is to study A-depcndcnce we 
must build our nucleus out of individual objects which do make sense for 
fixed coupling, and there are small dipoles. 

In the vicinity of the saturation boundary the amplitude for scattering of 
a dipole Q on a dipole p is given by [16] 

r^(/x, 6, Q) = T^{hp) ° InQ^IQl (8) 

where b is impact parameter and 

Ql{b,p,Y) = a{bp)p'[£n-a']—o " ■ (9) 
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To and a are of order one when b/j, is not large. Eq.(8) is valid in the region 



1 « inQ'/Ql « ^^x"(Ao)l^ (10) 

A-jr 

"^» 16(l-A„)3iV.,»(Ao) ^">^> (1^) 

where Aq is defined by x'(Ao) = -f^, and x(A) = ^(1) - |^(A) - i^(l - A) 
with ^0 the logarithmic derivative of the F-function. For simplicity of notation 
we now use Qs to represent the quark saturation momentum. Since the issue 
in this and the next section is the functional dependence of Qs the distinction 
between Qg and Qs is not important. 

The conventional way of looking at (8) is to view the dipole Q as mea- 
suring the gluon distribution of dipole /i so that 

dxG^{x,b,Q^) l^2fQ^^'~^' 



d?h a 



g^(|i) \nQ^/Ql (13) 



dxGu 



when Q'^/Ql > 1. (We note that when QVQs 

is on the order of one, ,->, - - 
Ql/a.) However, one may view the QCD-BFKL evolution in the opposite 
sense so that one considers the dipole Q to evolve to lower momentum 
scales, and then (8) gives the probability of finding a gluon (or dipole) at 
scale /X in the parent dipole Q. More precisely 

d .-..^ ..... 1 //^'A'"'\„./S2/ 2^ 



^^-7v^(g,6,y,/.)~-^c^j MQV/^') (i4) 

is the number of gluons per unit rapidity at scale ^ in the wavefunction of a 
dipole of size Aa;_L = 1/Q. In the large Nc limit would also represent the 
number of dipoles of size Ax^ > 1/// in a parent dipole of size 1/Q. In (14) 
we have introduced Qs which is defined by 

(161 

That is if one writes 
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Ql{b,fi,Y)=fi'f{b,Y) 

then 



(16) 



g^(6,g,y) = g7/(6,y), (i?) 

and 

T, r. a . (18) 

In Fig.l, the two directions of evolution are pictured. In the left-hand part 
of the picture the shaded saturation region for a dipole of size /i and rapidity 
Y is exhibited. Only the part of the saturation region where lnk\/ jj?' > 
is shown. In the right-hand part of the picture the part of the saturation 
region having lnQ'^/k\ < is shown for a dipole of size Q and rapidity Y. 
For the /i-dipole, evolution proceeds from y = up to y = Y, while for the 
Q— dipole evolution goes from y = Y — y = Qioy = Y. 

Now consider the scattering of a dipole Q on a nucleus made up of dipoles 
of size ji and nuclear density p. We view the process in the rest frame of the 
nucleus and in two steps. In the first step one takes the number density, 
of dipoles of size K and having rapidity y in the wavefunction of dipole Q, 
and in the second step one takes the scattering amplitude for a dipole of size 
K on the nucleus. Thus 



TA„ Q) . [ ,y 1^ . tAK, ,)f (19) 

where 

dN,{Q,y,K) ^ r dN,{Q,h',y,K) 
dy J dy ' 

The rapidity of the dipole K in (19) should be greater than InA^/^ so that the 
dipole is coherent over the size of the nucleus, but we suppose that InA^^^ 
is small enough that it can be neglected in setting the upper limit of the 
J/— integral in (19). Now 
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tA{K, b, n) ~ a\l/K^)£n K^fi^ ■ p2VR^ - (21) 

so long as tA is small, where we have taken the gluon distribution of a dipole 
H in the nucleus to be a in K"^ / fJi^ at scale K. Using (14), (20), and (21) in 
(19) one finds that the K integration (19) diverges in the infrared. But, this 
integration should be cut off at the value, Kq, when t^, given in (21) reaches 
one[21-24]. Thus 

TA-^Mr^'^n{Ql/Kl) (22) 

where, from (21), 

Kl ^ Q|(MTO ^ o?^n{Kl/^^)p2^W^^ I p\ (23) 



or 



^ ~ o?cA^lHn{p?cA'l'-) (24) 
A* 



with 
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{Kq is the quark saturation momentum in the Mclerran-Venugopalan model 
of a nucleus made of dipoles of size l/fi.) We shall always assume that 
a'^cA^/^ » 1 in order to have non-trivial nuclear effects. 
Now use (24) to eliminate Kq in (22). This gives 



Qi 

or 

The saturation momentum for the nucleus is defined as the value of Q"^ at 
which Ta becomes of order one. This gives 



//^ ^'^ ^■^ 1 — Aq 

where Ql{MV) is the saturation momentum in the McLerran-Venugopalan 
model. We note however that Ta does not quite have the usual scaling form 
(8) since 

Eq.(27) is only vahd for Q'^/Ql{A) > 1. When QVQK^) < 1 there is more 
than one dipole of scale Kq in the parent dipolc Q and unitarity effects will 
further suppress the Ta given in (27) [22-24]. One can express the scaling 
behavior of Ta more clearly by defining ariable 

Ql{A) = Q^^ca^A'/Hn{ca^A'/^) = QIQ^^{MV) / fj,^ (28) 
in terms of which, from (25), 

TA-^,[^y°in{QyQliA)) (29) 
where, now, (29) can be used when 
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>fi^4r. (30) 



Ql{A) V«' «'/ 
While Qs{A) gives the simpler looking scaling behavior it is Qs{A) which is 
the actual saturation momentum of the nuclear light-cone wavefunction. 
Finally, we evaluate Ta/T^j_. Prom (8) and (29) one finds 



1 (Ql{A)\'-'" IniQ^QKA)) 



Using (28) gives 

I-^/^i ~ ^j;^MV3)^*"'-^"(ca^^"3) ^1 _ _}__ j . (32, 

Thus not too far from the saturation boundary there is significant shadowing 
with the A-dependence of proportional to (A^/^)^^^". It is this shadowing 
which, according to Ref.ll, causes particle production in heavy ion collisions 
to scale,roughly, like Npart- We note that far from the scaling region, when 
^"q/^ >> 1) BFKL dynamics will replace Aq by 0, and an A^/^ behavior 
will again emerge for Ta- This is the region where double leading logarithmic 
behavior is valid. 



4 Running coupling BFKL dynamics 

In this section we revisit our discussion in the last section, but now using 
running couphng BFKL dynamics. We shall see that saturation looks quite 

different when running couplings effects are present. 

We begin by considering the scattering of a dipole Q on a dipole where, 
as before, Q/ ii S> 1. As in the last section[16-17] 



T^(/i,6,g)=.To(^ 



l-Ao 



1-Ao 



To is still a function of but now Qs is given by 



(33) 



Tih{l - Ao) 



Y 
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where 



„ jV,(l-Ao)|x"(Ao)]^ 

"=V — swAj — 

and 

c='-^ (36) 

with b = "~Y2^~~^5 and where A is the usual QCD parameter. ,^1 is the first 
zero of the Airy function, Ai[^). Of course, the form given in (33) is vahd 
only when InQ'^ /Ql 3> 1 so that the constant terms in (33) and (34) are 
arbitrary, and unimportant. 

The remarkable thing about (34) is that there is no yU— dependence present [16]. 
As noted earlier the /x— dependent corrections to (34) can be of size ^" ^^"^ ^ 
so that in case yU^/A^ ^ 1 there is a transition region between a low-F re- 
gion where fixed coupling dynamics occurs and a high-F regime where run- 
ning coupling dynamics occurs. This transition value of Y is[16] Ytrans — 

2hN~x{\)) '^{Ji)- ("^^) ^PP^y ^^^^ above this transition regime. 

The lack of a /i— dependence in Qa at large Y indicates the insensitivity of 
Qs to the nature of the target probed by the dipolc Q. Thus, (33)- (36) apply 
equally to a proton as to a dipole. And, of course, these equations must apply 
also to nuclei indicating that Qg, at large Y, has no A— dependence. This is 
in striking contrast to 74^/^-dependence of Ql found both in the McLerran- 
Venugopalan model and in fixed coupling BFKL evohition as given by (6) 
and (26), respectively. It is the purpose of the present section to try and 
explain why there is no A— dependence in Ql in running coupling evolution. 

It is useful to consider an expression for InQl/h? which interpolates be- 
tween fixed coupling evolution and running coupling evolution. To that end 
consider 



When r/£n2(/xVA2) < 1 
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p.^MA.VA') + '°';'rf'°' y (38) 
7r(l - Aoj 



while when Y / ln^{ij.'^ / K^) > 1 



- 7r6(l - Ao) ^ 2 / ^^-^(^°) y 

Y 7r6(l-Ao) 

SO that Ps matches onto the dominant parts of ln{Ql/K^) in both the fixed 
couphng regime, (38), and the running couphng regime, (39). Now view 

the cvohition running backwards, starting at a large scale Q and ending up 
at a smaller scale. As in the previous section call Qs the boundary of the 
saturation region in the backward evolution. Then in the approximation 
(39), and with ps — InQl/h?, one has 



p^ = 2ji^M^F 4iV.x(Ao) 

where p — £nQ'^/A'^. ps is, of course, a function of Q'^ and Y. Eq.(40) gives a 
property of the light-cone wavefunction of a dipole Q and may be applied to 
the scattering on any target, at least so long as ps remains in the perturbative 
regime. 

Begin by applying (40) to a large nucleus. Then, from (23), the scattering 
amplitude is at the edge of the unitarity limit (the edge of the saturation 
region) when 

where now replaces p^ in the running coupling case and for a realistic 
nucleus. Prom (40) 



determines Q^A). Now consider the scattering on a proton. Strictly speaking 
(40) does not apply since the proton is in the nonperturbtive regime. But, it 
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is clear that the unitarity hmit, for a fixed impact parameter, will be reached 
when Qs — A, that is when dipoles of size 1/A appear with high probbility 
in the wavef unction of the dipole Q. Thus 



determines Ql for the proton. Comparing (42) and (43) 



Y 7r6(l-Ao) 

where again we emphasize that QgiMV^ is the quark saturation momentum 
as determined in the McLcrran-Vcnugopalan model. Thus at large Y we see 
that there is no distinction between Qsi^A) and QgiP) although at current 
energies and for large nuclei the fixed coupling regime may be appropriate. 
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